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Abstract. A simple mechanims for the generation of elec- 
tromagnetic Zebra pattern emission is proposed. The mecha- 
nism is based on the generation of an ion-ring distribution 
in a magnetic mirror geometry in the presence of a prop- 
erly directed field-aligned electric potential field. The ion- 
cyclotron maser then generates a number of electromagnetic 
ion-cyclotron harmonics which modulate the electron maser 
emission. The mechanism is capable of switching the emis- 
sion on and off or amplifying it quasi-periodically which is a 
main feature of the observations. 

Keywords. Electron cyclotron maser, electron holes, auroral 
acceleration, auroral radiation fine structure, auroral kilomet- 
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1 Introduction 

Substantial effort has been invested into the explanation of an 
apparently miniscule effect in radio emission from some ob- 
jects, in particular from the Sun during type IV radio bursts. 
These radio bursts are caused in the aftermath of solar flares 
by electron-cyclotron maser emission from the anisotropic 
magnetically trapped flare-electron population. 

'Zebra' emissions are spectral fine structures which ap- 
pear as narrow-band multi-harmonic modulations (emis- 
sion/absorption structures) in the type IV radiation (for a re- 
cent review cf., e.g., [Chernov 2010| l. Their quasi-harmonic 
frequency spacing has sometimes been suspected to be re- 
lated to the ion cyclotron frequency, though in the past a 
number of mechanism have been proposed which search for 
explaining these fine structures differently (as for an example 
see, e.g.. 
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properties of the electron radiation or other processes like 
wave-wave interactions. 

Here we exploit the very simple idea that these emis- 
sions could indeed be a consequence of the presence of 
ion-cyclotron harmonics which are genuinely caused by an 
high energy ion distribution which propagates on the ra- 
diating electron background, itself radiating ion cyclotron 
waves which propagate on the R-X mode branch and, be- 
cause of their frequency w ^ lujd <^ uJce (with harmonic 
number £ = 1,2,3...) being much less than the electron cy- 
clotron frequency ujce, cannot leave the plasma. Interaction 
with the type IV electron background under certain restric- 
tive conditions then modulates the electron maser and causes 
the observed 'Zebra' fine structure. 

Justification of such a model (as is shown in Figure |4]i can 
be given by reference to the flare reconnection model, ac- 
cording to which flares are caused by magnetic reconnection 
in the lower corona at field strengths of several Gauss < i? < 
few 100 Gauss. In this process localised magnetic field- 
aligned electric potential drops develop which accelerate and 
heat the flare electrons and accelerate ions until they evolve 
into high energy essentially cold conical phase space distri- 
butions, a process well-knowm from auroral physics. Under 
type IV conditions the ions become relativistic. Their cold 
conical distributions excite ion-cyclotron harmonic waves 
via the ion-cyclotron maser instabiUty. 



2 Ion-cyclotron maser instability 

The cyclotron maser is a relativistic instability that involves 
the presence of at least weakly relativistic particles. It had 
been originally proposed for electrons only, first in the con- 
struction of free-electron lasers ( |Motz[ |1951[ |Madey| |1971| l 
and masers before being recognised for plasmas (by |Twiss| 
|1958| in its non-relativistic non-efficient version) and, much 
later (in its relativistic version by |Wu and Lee| |1979[l, as a 
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Fig. 1. Ion-cyclotron maser dispersion relation. Harmonic frequen- 
cies (blue) and growth rates (red) excited by a cold relativistic pro- 
ton ring distribution in an electron-proton plasma with ring-Lorentz 

(re-scaled data taken from Amato 



factor 7ii, < 10,N ii,/N ~ 10 

[and Aroris||2006^ . The ring-ion maser generates a large number of 
ion cyclotron harmonics £, all of nearly same maximum growth rate. 

very efficient radiation mechanism in space plasmas as well 



(for a review cf. Treumann 2006). Subsequently it has been 
identified as the main radio wave emission process in the au- 
rora of Earth and other planets, in solar radio bursts, particu- 
larly type IV bursts, and under astrophysical conditions. 

The electron cyclotron maser as the main emission mech- 
anism in all these objects is well established by now. Below 
we will make use of it. However, as for the Zebra fine struc- 
ture which interests us here, the electron cyclotron maser is 
just the carrier of information. Adopting the notion that the 
Zebra pattern maps the ion cyclotron frequency into the emit- 
ted electron maser radiation, we propose that an exactly sim- 
ilar mechanism works as well for the ions causing an ion- 
cyclotron maser instability. Such an instability has, in fact. 



been proposed first by Hoshino and Arons ( 1991 1 and was re 



fined by Amato and Arons | ( )2006j ) to include a cold relativis- 
tic gyrating ion beam. In a different non-relativistic version 



it was used in space plasma (cf., e.g., Chaston et al. 2002 1 
where it should serve to excite harmonics of low frequency 
electromagnetic ion cyclotron waves in the auroral plasma. 

For efficient generation of ion cyclotron harmonic waves 
under type IV burst radiation one requires the presence of a 
relativistic gyrating ion distribution with anisotropic velocity 
and positive perpendicular phase space gradient. For simplic- 
ity, the ion distribution is taken to be cold and anisotropic 
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where u is the normalised 4-velocity [u — 7/3 with 7 = (1 — 
/3^)~2 — (1 + ^2-^ 2 jjjg Lorentz factor], and Uu, is the initial 
value of u for the relativistic ions, and subscripts _L,|| refer 
to the direction of the magnetic field. For electromagnetic 
waves with perpendicular wave vector k = fcx the dispersion 
relation becomes 



-yx I ^x 



^xy ■ 
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= k^c^/u? (2) 



The explicit form of the dielectric tensor components is stan- 



dard and has been given, e.g., by Amato and Arons ( |2006| l 
whose results we will make use of here and who use a 
pair plasma as background. In order to exploit their results 
we have to re-scale their parameters to conditions of an 
electron-ion background. These parameters are the ratios a = 
B"^ / H^mNc^^ of magnetic to kinetic energy densities for the 
different species (in their case, the ion beam and background 
electron and positron pairs, index ±, for which they use 
numerically an, ~ 10^^,cr± = 2). We replace the positrons 
with ions and use the quasi-neutrality condition N = Ne = 
Ni + Nib, where the subscript ih denotes the relativistic 
ring ions. This yield for the background ion Lorentz fac- 
tor 7i = /Lt7e/ (1 — Nib/N) w /X7e. For small beam densities 
Nib ^ N, with /I — me/mi, the background ions are non- 
relativistic, and only the radiating electrons will be consid- 
ered relativistic. Moreover, we have ^ibhe ^ "^OOiJ^iN/Nib). 
Since the electron background is assumed mildly relativis- 
tic we will have jib ~ Q.1{N /Nib)je ^ 10 for applying the 
results of |Amato and Arons] (^2006 ). With these quite reason- 
able numbers we can directly refer to Figure 13(a) of Amato 
[and Arons] ( |20061 ) the re-scaled version of which which we 
reproduce here as Figure JT] 

One observes that a very large number of electromagnetic 
ion cyclotron harmonics is generated which propagate on the 
R-X mode branch. These waves have all roughly the same 
growth rate and are of narrow bandwidth. These harmonics 
are practically all confined to the plasma even in an under- 
dense plasma under the condition w^/cj^g < 2Tf.lm^(? < 1 



(Winglee 1983 1 of the electron-cyclotron maser instability 
as long as their frequency is below the upper hybrid fre- 
quency Lduh ^ <^ce in this case. The perpendicular electric 
fields E±{t) = ^££'£exp(— i^Wcii) of these ion cyclotron 
harmonics generate drift motions in the magnetically trapped 
anisotropic electron background plasma that is responsible 
for the type IV radiation. These electron drifts modulate the 
emitted electron cyclotron maser radiation, adding factors 
2TTS{ujzt£ujci) to the emission spectrum. 



3 Electron-cyclotron maser radiation 

We briefly indicate what kind of variation in the spectrum of 
the electron cyclotron emission is to be expected in this case. 
The model of type IV radiation is based on the assumption of 
a magnetically trapped electron distribution. Such distribu- 
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Fig. 2. Electron loss-cone (Dory-Guest-Harris) distribution in the 
normalised 4-velocity plane. Shown are two resonance circles for 
maximum and minimum radii, corresponding to minimum and 
maximum unstable electron cyclotron maser emissions in perpen- 
dicular propagation of R-X- modes. The section of the resonance 
circles contributing to emission are shown in red. Emission is pro- 
vided as long as the phase space gradient over the positive (yellow) 
part is steep enough to compensate for the larger part of the cir- 
cles traversing the blue negative gradient region of the distribution. 
Note that the blue region between the circles of maximum and min- 
imum resonances -Rmax,i?min absorbs electron-cyclotron waves by 
damping them. This heats the electrons in this domain and flattens 
the distribution in this region of phase space. 

tions are of loss-cone type which can, for instance, be mod- 
elled as (relativistic) Dory-Guest-Harris distributions 
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where v1 is the normalised thermal spread, and the factor 
v?f oc sin^'* 9, with 6 the angle between 4-velocity and mag- 
netic field direction, accounts for the loss cone. Note that s 
can be any real number while being most conveniently cho- 
sen as s G N. 

With the above distribution function the instability has 



been exhaustively discussed by |Pritchett| ( [1986) for differ- 
ent propagation angles finding that largest growth is obtained 
for strictly perpendicular propagation, fc|| = 0. In this case 
the resonance condition for resonant excitation of waves be- 
comes Ti'^ + ?i| — 2(1 — i^ce) = 0, where i^ce = ^/oJce is the 
non-relativistic ratio of wave frequency to electron cyclotron 
frequency. Clearly, resonance is possible only for i/ce ^ 1 
and, as usual, the resonance Une is a circle in the normalised 
4-velocity plane {u±,un ) of radius 



(4) 



located between the minimum and maximum radii 
Rmin,Rmax of resonance (see Figure [2| which corre- 
spond to maximum and minimum resonant frequencies, 
respectively. 

The relativistic cold plasma dispersion relation of the R-X 
mode, which is the real part of 



l-(2-ni)A = 0, with 



allows for a range of such resonant frequencies below the 
non-relativistic ujce which is quite narrow (cf., e.g., Pritch- 
ett 1986 1, depending on the plasma parameters. The max- 
imum resonance frequency (minimum radius) is very close 
to Vce = 1- On the other hand, the minimum resonance fre- 
quency can be estimated from the cold X-mode dispersion 
relation yielding 
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The function A contains the plasma response. Taking fc|| = 
0,|1 — i^cel ^ 1. one has to first order Ke{A) ^ 1. This ex- 
pression is to be used in the calculation of the growth rate 



Im(j/ce)=i-ilm(^) 



(7) 



Maser radiation will be produced if lm{A) < is negative 
(corresponding to 'negative absorption' of the electromag- 
netic R-X mode). The resonant imaginary part of A is given 
by 



Im(yl) : 
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In dealing with the loss cone distribution Figure [2]one re- 
alises that the modulation of the particle distribution caused 
by the presence of the ion cyclotron waves just affects the 
low perpendicular energy electrons at the loss cone bound- 
ary. This is the only region which has positive perpendic- 
ular gradient and contributes to wave growth over the nar- 
row part of the boundary crossed by the resonant circle. The 
longer part inside the distribution sees a weak negative gradi- 
ent and causes wave absorption which for instability is over- 
compensated by the steep gradient at the loss-cone boundary. 

Inspection of this figure suggests, that the modulation of 
the loss-cone boundary will heavily affect the emission con- 
ditions. In presence of ion-cyclotron waves the low perpen- 
dicular energy electrons at the loss-cone boundary will start 
oscillating up and down in u± with the imposed frequency 
of the ion cyclotron wave electric field, when performing an 
Ej^ X B-drift. This drift displaces the loss-cone boundary ei- 
ther up to higher u± or down. In the latter case the loss-cone 
is partially filled, and the maser may switch off, causing zero 
emission intensity at the respective frequency. In the former 
case the loss-cone widens, faking a larger loss-cone exponent 
s and bringing the minimum-radius resonance almost com- 
pletely into the positive gradient region which should am- 
plify the emission at frequency very close to i/ce ^ 1. These 
two cases are shown schematically in Figure [3] 

Including this modulation effect into the calculation of the 
electron cyclotron maser instability from scratch on is quite 
involved. We can however boldly construct a simple model 
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Fig. 3. Same as in Figure |2] except for tiie modulation of the loss 
cone included. Top: The modulation due to the electric field drift in 
the transverse electric field of the ion-cyclotron mode shifts the per- 
pendicular velocity up thereby opening the loss cone and pushing its 
gradient to higher u±. This steepens the gradient and reduces the 
damping part of the resonance circle. At minimum radius (largest 
emission frequency Uec) the resonance circle is fully in the positive 
gradient. During this phase emission is amplified. Bottom: During 
this phase the drift reduces the perpendicular velocity filling the loss 
cone and switching the maser off. 

which possesses some of the relevant features of the modu- 
lation. 

Realising that the two time scales of the ion cyclotron and 
electron instabilities are quite far apart, we consider the mod- 
ulation as stationary on the time scale of the electron cy- 
clotron maser emission. Then, the drift effect will result in 
a modulation of the loss cone as has been described above 
and can be included by replacing u± -^u'j_=u± + A in the 
loss-cone factor of the electron distribution function Eq. Q, 
where A = jeVE{t)/c is the contribution of the normalised 
drift in the ion-cyclotron wave electric field. 

Modulation can, in addition, also be included as a varia- 
tion of the loss-cone exponent s replacing it with a variable 
exponent s -> a{s,t) > and t the slow time scale, allowing 
for a direct response of the loss cone to the modulation. Since 
a is not anymore an integer number, this latter step requires 
that in the normalisation of the DGH distribution Eq. ^ the 
factor si is replaced by the Gamma-function r(a + 1). Per- 
forming the derivative yields 



duj_ 



(9) 



to u^, rearranging and neglecting terms higher than first or- 



Inserting into Eq. ([8]), performing the integration with respect 

than f 

/ a 



der in A reduces to 

Im(yl)«- 



r(a+l) i/ceW^g w2 



x°'dx 



A' 

R'res 



x-i +A'(2a+1) 



(10) 



where A' = + R^es (t) / c. All quantities are nor- 
malised as indicated earlier. (In order to arrive at the above 
expression we have taken Rres out and changed the fake 
integration variable.) The integral can be evaluated mak- 
ing use of J^dx'^/\/l — x = B(a + 1,|), where B{x,y) = 
T{x)T{y)/T{x + y) is the beta-function and T{^) = y/n. 
This yields 



Im(A)«-C(0i+A'Q2) 



(11) 



where C = {7rRresU}^)e''^^" /(I'cei^^e'^l) > 0. The term Qi 
depends solely on the loss-cone exponent a, while Q2 de- 
pends on a and accounts for the variation in the perpendicular 
velocity. These functions are 



Qi = 



1 



-1 



(a+i)r(a+i) 

a(2a + l)r(a+5) 
(a + l)r2(a + l) 



i?3 
"res 



Rr 



Observing that always a > we first neglect Qi. Then, Q2 
[and Im(yl)] changes sign at a = R^J-^ = [2(1 — Vce)]"^ ■ If 
a = s is an integer, this means that at sufficiently large loss 
cones and for A' > the term with Q2 damps the resonance 
and could theoretically cause absorption. On the other hand, 
since |;/ce — 1| <C 1 this requires quite large a. We reasonably 
assume that a < R^^^ and thus Q2 > 0. However, the sign of 
this term also depends on A' oc oc sm{£ujcit) which os- 
cillates with harmonic ion cyclotron frequency £ujci- As we 
have proposed, this causes a sinusoidal modulation of the 
term including Q2 and leads to modulation of the electron 
maser emission at the resonance frequency. At small a, how- 
ever, we have always Q2 <Qi- In order for the second term 
in Im(yl) to dominate one requires that with /3e = | V^^^/cj 

(2ai?L7e/3E)"' < ^ ~ |sin(£w„t)| < 1 (12) 



where, for simphcity, we assumed just one sinusoidal ion- 
cyclotron harmonic belonging to the entire spectrum of ion 
waves that are generated by the cold ion conic distribution. 
The argument of siii^^ must be smaller than one. Hence this 
condition, combined with the condition on a 



Rres>a:>{2jel3E) ^ 

and with 7e > 1 implies that 

Rres < We or l-Uce<. 2/9| 



(13) 



(14) 
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Fig. 4. The model of a post-flare magnetically trapped electron dis- 
tribution which is the source of type IV electron cyclotron maser 
radiation. In this stage the weakly relativistic hot electrons have 
evolved into a loss-cone distribution of Dory-Guest-Harris type 
with completely empty sharply edged loss cone assuming that 
whistler activity has ceased (simply because ujpe/^cc ^1) and 
quasilinear diffusion into the loss cone is terminated. The R-X mode 
emission is caused on the section of the resonance circle in phase 
space that cuts through the steep perpendicular phase space gradient 
of the loss cone and is mostly perpendicular to the magnetic field 
at frequency below but very close to the local electron cyclotron 
frequency. Upward directed strong electric fields at the site of the 
magnetic mirror accelerate cold ions into an ion-conic distribu- 
tion which propagates up into the configuration and generates large 
numbers of ion cyclotron harmonics. The local modulation of the 
electron loss cone distribution by the ion wave spectrum modulates 
the electron cyclotron maser and causes the genuine ion-harmonic 
spaced Zebra emissions. 



These conditions can be satisfied only for sufficiently steep 
loss cone gradients a and large drifts caused by the ion- 
cyclotron harmonics. If this is the case, the growth rate Eq. 



(Ill of the electron-cyclotron maser instability is periodi- 
cally modulated during the phases of the ion-cyclotron har- 
monics. It is amplified during positive phases < uidt < it 
and switched off during negative phases tt < ujcit < The 
frequency of this modulation of the electron-cyclotron maser 
radiation at the local resonance is proportional to the ion cy- 
clotron frequency lOd- 

Though the modulation is small in the growth rate, it mod- 
ulates the exponential amplification factor of the maser emis- 
sion which causes a much stronger modulation of the emitted 
radiation intensity. However, it requires a sufficiently rela- 
tivistic trapped electron component with steep loss cone gra- 
dient in addition to the presence of a cold relativistic ion ring 
passing across the electron background. 

In order to directly see the modulation, we calculate the 
power emitted by the electron-cyclotron maser at resonant 



1 



(15) 



where we have assumed that utcei^ce/k ~ c. 

We are interested only in the modulation of the power We 
therefore do not calculate the power at saturation which re- 
quires the precise knowledge of the saturation process. Since 
in the collisionless case the radiation can presumably freely 
escape this requires knowledge of the competition between 
cyclotron maser emission and cyclotron re-absorption along 
the ray path. This is not impossible to estimate but lies out- 
side the scope of our investigation. So we stay with the linear 
state, assuming that the power is limited just by escape from 
the volume, i.e. by the time t^sc ^ L/c, where L is of the 
order of the linear dimension of the trapped electron volume 
perpendicular to the magnetic field (see Figure |4]|. 

The first term in the exponential accounts for the am- 
plitude the power can reach, the second term accounts for 
the cyclotron modulation via the time dependence of A' cx 
sm{£Ldcit), where the time t is now measured on the scale 
of the ion-cyclotron period. Evaluating the exponential term, 
we find for 



,,=1 



£uJrd 



(16) 



where /^(z) is the modified Bessel function of order rj, 

and z = 2TescCQ2(3E\/^ + Rres - Oils observes that the term 
Io{z) on the left does not produce any modulation. Since the 
power is a real quantity, only the real part of this expression 
contributes. Thus we obtain 

^(i^ce,re,e,i)=7'o|/o(^) + £/r,(z)sin[r;(fa;,,i)]| (17) 
where (€,77) S N, and we have written 



Vo = \E^{iy,e)\l exp[2CQiTesc] 



(18) 



These expressions show indeed that the emitted maser power 
experiences a modulation at multiples of the ion cyclotron 
frequency. In frequency space ro, defining ^ — iujci, this cor- 
responds to 

V{vce,i^) oc 2ttVo{Io{z)5{zu)+ 

+ J2\I,{z)[S{w-TjO-S{w + vm \ (19) 

r,= l J 

The first term on the right describes just the stationary spec- 
trum while the sum superimposes the ion cyclotron modula- 
tion on the spectrum. 

The above mechanism, a sketch of which is given for il- 
lustration in Figure |4] may occasionally work and produce 
the wanted Zebra structures with periodic amplification and 
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suppression of electron cyclotron maser emission in type IV 
radio bursts. It requires a rather steep or abrupt transition be- 
tween the electron distribution and the loss cone. It, how- 
ever, also depends on the excitation of strong ion cyclotron 
harmonic waves in order to satisfy the restrictions imposed. 
This assumes that an intense ion (conic) ring distribution has 
been generated during the flare, a process that works only 
in the presence of sufficiently strong field aligned electric 
fields produced in a magnetic mirror configuration. Because 
of these reasons it might not be realised very frequently. 
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